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I. INTRODUCTION 



Superpartners of gauge and Higgs fields play an important role in SUSY phenomenology. 
In particular, neutralino dark matter (DM) in the SUSY framework was considered in detail 
(see, for example, Refs. and these investigations were adapted to astrophysics. So, 

an analysis of the neutralino system and the structure of gauge bosons interaction with the 
neutralino and chargino is important for the DM description and the study of astrophysical 
data. 

In many phenomenological works both the neutralino mass spectrum and the structure of 
states follow from the formal diagonalization of the neutralino mass form by an orthogonal 
(real) matrix [H, [sl - Such a procedure does not consider some important features of 
the structure of the Majorana states, related with the sign of the mass. These features are 
connected with the structure of the neutralino-boson interactions which, in turn, defines the 
peculiarity of the neutralino-nucleon scattering. 

The most cornplete and comprehensive analysis of the neutralino system has been per- 
formed in [12I - 3]. In these papers, the diagonalization of the neutralino mass matrix is 



considered in detail in the MSSM and some of its extensions. Special attention was paid to 
the building of neutralino states with positive masses. However, due to the complexity of 
the general diagonalization formalism it is difficult to trace a link between the sign of mass 
and the structure of the neutralino-nucleon interactions. 

In this paper, we analyze the features of the neutralino structure and interactions which 
are directly related with the sign of mass. We consider the simplest case when this con- 
nection is transparent and convenient for illustration. In the second section we compare 
two ways of the diagonalization - by orthogonal and unitary diagonalyzing matrices. These 
two variants lead to neutralino states with opposite and equal signs of masses. They are 
formally equivalent and related by a field redefinition (see Section 2). But the negative 
mass of the neutralino (as it occurs for the first case), has to be taken into consideration 
in the consequent calculations. Disregarding this important feature, it is possible to get an 
incorrect conclusion on the spin-dependent (SD) and spin-independent (SI) contributions 
into the neutralino- nuclear cross-section [lo| . 

Redefinition of the field reveals a link between the sign of the mass and transformation 
properties of Majorana spinors with respect to inversion (i.e. parity). An analogous connec- 
tion between the sign of the mass and parity was revealed for the case of massive Majorana 



neutrino in Refs. [171] and [18l]. In Section 3, we present the compact Lagrangian of the 
boson-neutralino-chargino interaction in terms of a redefined field, which is convenient for 
phenomenological applications. 

In the fourth section, we consider the neutralino mass matrix diagonalization by means 
of a unitary matrix giving all positive masses. Thus, the standard calculation rules can be 
kept unchanged, and there is no need to check the sign of the mass or redefine the field. A 
convenient diagonalization procedure based on the exponential parametrization of the uni- 
tary matrix is discussed. This procedure is formalized in a perturbative calculation scheme 
analogous to 13] ■ However, our scheme needs a smaller number of input parameters and 



gives all expressions in a quite compact form, which is useful for calculations. The method 
suggested is generalized for the case of mass matrix with complex parameters (Appendix 
B). 
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II. NEUTRALINO PARITY AND STRUCTURE OF BOSON-NEUTRALINO IN- 
TERACTION 



Now we analyze the connection between the signs of the neutrahno mass and the struc- 
tures of the neutrahno-bozon interaction when Mz/M^. — > and is Mi, M2 or /x. In this 
hmit, the analysis is simplified considerably, but the results can be used in the general case 
too. This limit is approximately realized in Split SUSY scenarios [sj and strictly takes place 
at high temperatures T 3> Eew-, when the Higgs condensate is melted (the high symmetry 
phase). For completeness, we give the well-known minimal formalism that we need in the 
following analysis. 

If the mixing of gauge and Higgs fermions is neglected, the mass term of higgsino-like 



Majorana fields has the Dirac form [19 
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Mk = ^l^mnm^ + mRHlL) + h.c. . (2.1) 

This form can be represented by a (2x2)- mass matrix which is known as the specific 
matrix with zero trace: 

M, = (° ^) . (2.2) 

There are two ways to diagonalize this matrix. The formal procedure using the orthogonal 
matrix O2 leads to a spectrum with opposite signs: 

O^M^O^ = ((J , O2 = ^ (} _\) , = (/X, -^^) , (2.3) 

where Tr{OjM202} = Tr{M2} = (trace conservation). In this case, one of the Majorana 
fields has a negative mass, regardless of the sign of /i. The matrix M2 can also be diagonalized 
by the unitary complex matrix U2, giving masses with the same sign: 

U^^^^U, = °) . U. = -L (1 ^^) . = (2.4) 

The diagonalization (12.41) is equivalent to the procedure (12. 3p with the redefinition x ~^ ^75X 
of the non-chiral (full) field with m = — /i. The last transformation is equivalent to Xr,l 
±iXR,L for the chiral components. Note that there is an infinite set of unitary matrices 



= U2 ■ which diagonalize the mass matrix M2 (see also [15|, Appendix A. 2): 



1 / e^'t' ie^'t' 



g.v' ^g'v \ _ / COS — sm 



~ ^2 (e~** -le-'^ J ' - (^sin0 cos /J ' ^^'^^ 

The additional O2 - symmetry (see Appendix B) leads to a free parameter arising in the 
general case. 

Dealing with the spinor field we should take into account the sign of its mass in the 
propagator and polarization matrix or redefine the field with a negative mass. As a rule this 
feature is not considered in phenomenological applications (see, for example, Refs. 11|). 



From the redefinition x' = ^TsX? it follows that the transformation (relative to inversion) 
properties of Majorana fields having opposite mass signs are different. As a result, we have 
one usual Majorana field and one pseudo-Majorana field. 
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The gaugino mass subform is of the standard Majorana type [19| and has no specific 
features. The signs of the masses for xi and X2 are defined by the signs of Mi and M2 
in the case of small mixing. They can be made positive by a redefinition. Note that the 
redefinition procedure always influences the mixing terms of the mass matrix and should be 
taken into account in the general case (see Sect. 4). 

Now we consider the connection between the structure of the boson-neutralino interaction 
and the relative sign of the neutralino masses. For simplicity, we show this connection in 
the pure higgsino approximation. The contribution of terms caused by mixing is considered 
in the next section. We here present a short comparative analysis of the calculation rules 
in two cases: when the masses of xs and X4 have different signs (diagonalization (12.31) ) and 
when they have the same signs (diagonalization (I2.4p ). The initial Lagrangian is 

hut = \gzZ,{m^^^Hl^ + ^2V^2r) , (2.6) 

where gz = 92/ cos 6w- The diagonalizations (12. 3p and (12.41) lead to the following forms of 
neutralino-boson interactions, respectively: 

1 / i 

(1) Lint = -^gzZ^X3l''l5x'i; (2) Lint = ^gzZ^.Xsl^X'i- (2.7) 

In Eqs. (12. 7p the first case with opposite signs (/i, —fi) can be transformed into the second 
case with the same signs (/x, /i) by the redefinition ^75X4 = Xi- Here we show that both 
Lagrangians in Eqs. (12.71) lead to the same result without any field redefinition if the neg- 
ative sign of X4 mass is considered in the calculations. In other words, both structures in 
Eqs. (12. 7p lead to the parity-conserving vector interaction which gives the spin-independent 
contribution to the neutralino-nucleon scattering poj. 

Let us consider, for example, the process of the scattering xsQ X'iQ with t-channel 
exchange of a Z-boson (t^ <^ M^). In both cases, the amplitudes of this process are 

(1) Ml ~ xtiP2h'"l5X3iPi) ■ tihhf^icg - l5)q~{ki) , 

(2) M2 ~ xUp2)i^Xz{Pi) ■ tihh^ic, - i,)q-{h) . (2.8) 

Formally, the amplitudes A4i and M.2 have a different structure. Therefore, one can drow 
a wrong conclusion about the contributions to the spin- dependent and spin-independent 
parts of the cross-section , if the negative sign of Xa mass has not been taken into account. 
However, taking into account the negative sign of Xa m the polarization matrix allows one 
to get the same result for both cases. If x is in an initial or final state, the polarization 
matrix of the field x in M.^M. is defined by (for positive mass = /i > 0) 

E^^(P)^-(P) = ^(^±/^)' (2.9) 

or (for negative mass = — /i) 

E Xl^{p)xl^{p) = ^{pTfi). (2.10) 

(T 

With the help of Eq. (^UO^ . we get 

MtMi ~ Tr{(p2 - /x)7'^75(Pi + = Tr{(p2 + /x)7^(Pi + (2.11) 
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One can get the same expression for A^^A^2 using the standard definition (12. 9p of the 
polarization matrix. This feature should be included in an analysis of neutralino-nucleon 
scattering. From the interaction Lagrangian only, without consideration of the mass signs, 
we cannot drow any valuable conclusions on the SD or SI contributions. In particular, 
the bilinear structures X3l^lXi Xsl^lnXi are vectors, while X31^lXi and X37^75X4 are 
axial vectors. Analogously, X3X4 and X375X4 are pseudoscalars, while XslbXi and X3X4 are 
scalars. Thus, the analysis of the neutralino-nucleon interaction has to take into account 
neutralino transformation properties. As a rule, in the bulk of papers this feature has not 
been considered explicitly and mistaken conclusions can be obtained in calculations of the 
SD and SI contribution to the neutralino-nucleon interaction. In particular, for the current 
structure Xil^l^XkZ^ it is possible to obtain SD or SI neutralino-nucleon cross sections 
depending on the neutralino relative parity. For instance, in Refs. |2l|, j22|] the same 
current structure was considered without any comments on this important peculiarity. From 
our analysis, it follows that in the case discussed, the neutralino-boson interaction gives the 
main contribution to the spin-independent part of the cross-section . 



An analogous feature is in order when Xi is in an intermediate state, for example, in the 
process X3^ ^ X4 ^ X?,Z- The amplitude Mi of the process is 

Ml ~ xt{P2)l^l^{q - /w)7"75X3"(PiXef = xt{V2)l^{q + m^)Y Xz iViy^u ■ (2-12) 

In Eq. fl2.12p we use the propagator ~ — /u) for the field xl with negative mass = — /i, 
whereas the standard propagator is ~ (g + /i). So, the mass sign been taking in account 
leads to the same result for the amplitudes TWi and M.2, where A^2 describes the same 
process with redefined X4 iu an intermediate state. 



III. GAUGE BOSON-NEUTRALINO-CHARGINO INTERACTIONS 

In this section, we give compact expressions for the Lagrangian of gauge bozon-neutralino- 
chargino interactions in the case of small mixing. These expressions are convenient for 
calculations in cosmology. This Lagrangian follows from Eqs. (I5.14p - fl5.16l) (Appendix A) 
as a result of the shift: 



Hut = '-g2eabcW''rW'W'^-^g,H{YH-^B,+ 

+ ^g,W;{H^^,Hl^ + HlYHti) + (3.1) 

+ ^g^Wli-H^^^H-^ + Ht^Htj^ + m^Hl^ - 

- mrH^2L) - lgiHh'H',,B, + ]^giHh'Hl^B,. 

Let us consider the case Mz -C /i. Mi 2, which can be used in Split SUSY models [2|, 0], 
0], 10]. The physical states of the neutralino in the zeroth order of the mixing were defined 
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in Section 2, and the chargino states in Appendix B: 



X4 = ^75(^^l° - i^2)/v^; (3.2) 



In Eqs. (13 ■2p we do not use charge sign notation for the Dirac fields H and W (in contrast 
to in analogy to the Standard Model notation. From the structure of H in (I3.2p . it 
follows that the components if^^ and B.^^ correspond to particle and anti-particle parts in 
a Weyl basis. Using the definitions (13.21) we represent Lint in the form 

+ 92W;{W^^xi + \h^'x^ - \hYxa) 

- g2 COS OwZ^Wj^W - ^ ^\ cos 2ewZ^HYH (3.3) 

2 cos Ow 

2cos6'vi/ 

The first order corrections to the ZxiXk interaction caused by the mixing (see Appendix 
B) are 

L'Jix =7> ^/^ XiTX?. + T7 X2TX^ 

2 cos Uw Ml — /i M2 — /i 

-Xi7^75X4 + 7^-^X27^75X4), (3.4) 



where the are defined by Eq. (14.51) in the next section. From Eq. (13.41) . one can see 
that the interactions of X3 and Xi with xi,2 have a different structure. This effect is directly 
connected with different signs of the masses of the non-redefined fields. Note also that 
the bino-like neutralino X2 ^ B does not interact with gauge bosons in the zero mixing 
approximation, but it interacts with the scalar Higgs field and Xsa- 

In the pure higgsino limit, Xs and X4 constitute the neutral Dirac field = (X3+^X4) / 
and the part of Eq. (13. 3p can be represented in the form (here we omit the heavy states W 
and X2) 

+ ir^zfi'^.H' - (3.5) 
2 cos Ow 

- ^ ^\ cos2ewZ^H^^H-eA^HYH. 
2cos6'i4/ 

The Dirac representation (13.51) of the boson-neutralino-chargino interactions involving a 
small mixing of the Higgs fermion with the gauge ones is formal (unphysical) but is conve- 
nient for our calculations. In this case, we avoid some complications of the Feynman rules, 
caused by the Majorana nature of xs and X4 23] - 2^. By direct calculation we have 



checked that both ways lead to the same results for the annihilation and co-annihilation 
cross-sections UL 
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IV. DIAGONALIZATION OF THE NEUTRALINO MASS MATRIX BY UNI- 
TARY MATRIX WITH EXPONENTIAL PARAMETRIZATION 



In this section, we consider diagonalization of the 4x4 mass matrix with real parameters 
fi, Ml, M2. Generahzation of the approach for a matrix with complex parameters is consid- 
ered in Appendix B. As follows from Section 2, the sign fi is not essential, and Mi 2 can be 
made positive by a redefinition of the gauge fermion. The neutral fermion mass form follows 
from the SUSY Lagrangian (Eqs. flS.lSp and fl5.16p in Appendix A) after the shift 



1 

■2' 



^R) 



M, 



o<Pl 



+ h.c, 



(4.1) 



where 



[B, H^, Hi) and 



Mn 



/ Ml 


-iMzSeCp 
\ iMzSeSp 



-iMzSeCp iMzSeSp \ 

M2 iMzCeCp -iMzCeSp 
iMzCeCp 

-iMzCeSp fi / 



(4.2) 



where sg = sin 6 and = cos (3. The matrix (14. 2 p differs from the commonly used one by 
the presence of the imaginary unit in the mixing terms. One can go to a real traditional 
matrix Mq by a redefinition if^ = i'j^Ha, where a = 1, 2. As a result we get the standard 
matrix following from (14. 2 p under the formal transition iMz Mz and fi —fi. However, 
impl ying our diagonalization procedure there is no need to do this transformation (see also 
Ref. H). 

It is convenient to analyze the diagonalization of the matrix (14. 2 p with the help of the 
intermediate transformation 



M, 



1 

U2 



(4.3) 



Here 1 and are the identity and zero (2x2)- matrices, and U2 is defined by Eq. (12.41) in 
the pure higgsino limit. Then the intermediate mass matrix has the form 



M, 



where 



mi 
m2 
ms 
m^ 



( Ml —imi m2 \ 
M2 imj, —m/i 
-imi im^ fi 
y m2 —m^ fi J 



Mz sin ^vy(cos/3-sin/3)/V2, 
Mz sin Ow^cos f3 + sin (3) /v^, 
Mz cos ^vi/(cos/?-sin /3)/V2, 
Mz cos 9 w {cos P + smP)/V2. 



(4.4) 



(4.5) 



Intermediate fields are defined by 0/ = that is, (0/)^ = {B,W^, X3,xi)^ where xi 

and xi ci-re defined by Eq. (3.2). The use of the intermediate mass matrix provides the 
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positivity of the higgsino-like neutralino masses and leads to the " quasidiagonal" structure 
of the matrix in the case of small mixing. 

The matrix fl4.4p is symmetric and complex, but it is not Hermitian. The spectrum of M/ 
is real and has a simple form. However, it is not a mass spectrum of the neutralino, because 
the diagonalization of the neutralino mass matrix UjM/U/ = diag{mf^) differs from the one 
defined by U+M/U = diag{Xk)- In the last case, U is built of the eigenvectors of M/ and 
the Xk are the eigenvalues of M/. According to the theorem 4.4.4 (Takagi expansion) from 
Ref. 27|, any complex symmetric matrix can be diagonalized by the unitary matrix U: 

U+MU* = diag{mk), > 0, (4.6) 

where U is built from eigenvectors of the matrix A = MM* with the spectrum {ml}, i.e. 
U^AU = diag{m\). Consistency of the last relation and Eq. (14.61) is evident from the 
equality 

U+MU*(U*)"^M*U = rfm^(m^), ml = mk. (4.7) 



The method based on the Takagi theorem was considered in Refs. [12j, [15j, [28j and |29| . 
where the standard way of the determination of the spectrum is given. However, there is 
no need to solve this complicated problem in the case considered. Here we show that the 
spectrum of the matrix A = M/M^!" coincides with the squared spectrum of the traditional 
real mass matrix Mq. The spectrum of the matrix A = M/M^ (M/ is defined by (14.41) ) 
follows from the solution of the characteristic equation det(A — A ■ 1) = 0, 

- aX^ + h\^ -c\ + d = Q. (4.8) 

The coefficients a, 6, c, and d in Eq. (14.81) are expressed in terms of the matrix elements of 
Mq as follows: 

a = Ml + Ml + 2ii'^ + 2Ml; 

b = MlMl + 2Ml{n^ + M| cos^ Ow) + 2M2(^2 + 

+ M| sin^ Ow) + 2M|/i sin 2/?(Mi sin^ Ow + 

+ M2Cos'ew) + {^^ + Mlf- 
c = 2n'^MlMl + M2[(/i2 + M| cos^ Owf + 

+ 2M2MI/X cos^ Ow sin 213] + Mj [(/x^ + M| sin^ Owf + 

+ 2MiM|/i sin^ Ow sin 2(3] + ^MiMaMf sin^ 2dw + (4.9) 

+ 2Mln^ sin 2/?(Mi sin^ Ow + M2 cos^ Ow) + 

+ ;U^M| sin^ 2/3; 

d = i/MjMl + M2/i^M|(2/iM2 + M| cos^ Ow sin 2/?) x 

X cos^ Ow sin 2/5 + MliJ,'^Ml{2fiMi + M| sin^ Ow sin 2(3) x 

X sin^ Ow sin 2(3 + -MiM2/i^M| sin^ 2ew sin^ 2(3. 



Analogous expressions are given in Ref. 13], where the algorithm of the definition of the 
spectrum Afc is considered. In the general case, we can get exact expressions for the roots A^ 
of Eq. (14.81) in terms of its general algebraic solutions. It is difficult to analyze and compare 
such expressions, but we can show that the roots of Eq. (14.81) are A^ = ml, where m^ is the 
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conventional neutralino spectrum. To show this, let us write the characteristic equation, 



det(MQ — / ■ 1) =0, in the form (Mq is the standard real mass matrix [291], 31 1) 



(Ml - /)(M2 - /)(/^ - /i^) + M|(/ + /isin2/?) x 
X {Micas'^ Ow + M2sm'^ew -I) =0. (4.10) 

Then we arrange the even and odd degrees of / on the left and right hand sides of this 
equation separately. Squaring the equation and changing ll A^, we get the equation 
(14. 8 p with coefficients (14.91) . Moreover, by direct calculation we have checked that the mass 
spectrum appearing as a result of the diagonalization 

V^MiV* = diaginik) (4.11) 

is entirely positive (see Appendix B). 

With the help of the Takagi theorem it is possible to illustrate the correct construction 
of the positive mass spectrum. However, the above discussed method is not convenient for 
calculations and can be used as the diagonalizability proof only: there is a unitary matrix 
with the property (14. 6 p or (14.111) . The use of Mj gives a convenient tool for the calculation 
of the spectrum and states when Mi,M2,fi and their differences are much greater than 
Mz- The hierarchy of Mi, M2 and /i is arbitrary, i.e. one can apply the method suggested 
to the various scenarios of the neutralino DM. In this case, the diagonalyzing matrix is 
quasidiagonal, i.e. \Ukk\ ~ 1 and \Uik\ <^ 1, i ^ k. Then we represent this matrix in the 
exponential form which contains six angles and six phases as input parameters. A similar 
approach was considered in the general case in Ref. [12], where six angle and ten phase 
parameters were applied. Here we show that in the case of the mass matrix (14.41) it is 
possible to use six phases only (see Appendix B): 



U 



rie*°i a2 5^6-''*'^ S^e'''!'^ 



(4.12) 



where Sk and (pk are input angle and phase parameters. The values a/3,rj and are some 
functions of the input parameters which are defined by the unitary conditions {\J~^\J)ik = 
Sik- In our case, \Sk\ -C 1 and functions ap,ri and ak are easily determined by successive 
approximations (Appendix B). Apparently, the diagonalization of the real matrix demands 
the angle parameters only. Having used the diagonalization conditions 

Uf M/U/ = Mrf = diag{mi, m2, m^, 7714) (4-13) 

the input parameters Si and (pk can be determined from the six independent equations 
(Mrf)ifc = 0, i > k. So, there are six conditions for the real and six ones for the imaginary 
parts of matrix elements. Then the masses rria appear as functions of the defined input 
parameters. As is shown in Appendix B, the perturbative calculation scheme can easily be 
formalized. 

The above discussed method of diagonalization is applied to the case of a mass matrix 
with complex parameters Mie*'^^, M2e*'^^ and /ie*'^'" (see Appendix B). In this case, we 
have to generalize (14.121) introducing additional phase parameters according to ak — > ake^^'". 
Thus, we have the same quantity of parameters as in [l2|. The functions ap^Vk and ak 
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are determined in terms of the input parameters Sk, 4>k and using the unitary condition 
U"'"U = 1. The input parameters are determined in terms of the mass matrix elements if 
we use the diagonahzation conditions (U^MU)ifc = 0, i ^ k and Im(\J'^'M\J)ii = 0. Note 
that in this case the perturbative calculation scheme (as for the real mass matrix) is in order 
also. However, the expressions are more complicated and bulky, so we give the results in 
the first approximation only (Appendix B). 

Here we represent the mass spectrum and parameters of the matrix U defined by (14.121) . 
up to terms ~ m|/M^, a = 1, 2 (Appendix B). The neutralino masses are 



^X2 



Mi + ^^|r^(M,+;.sin2/3), 

IVlo /i 



M|l-sin2/3 

^ 2(Mi + /i)(M2 + ^) ^ ^ 

X (Micos^6'w/ + M2sin^6'w/ + 
M|(l + sin2/5) 

X (Ml cos^ Ow + M2 sin^ Ow - f^)- 

From Eqs. (4.14) one can see that 777.3 aiid have the same sign. The validity of the 
expressions (14.141) does not depend on the hierarchy of Ma and fi. So one can use Eqs. (I4.14p 
in various SUSY scenarios. 

Another feature of the diagonahzation is the presence of a free parameter in the structure 
of neutralino states (Appendix B). Evidently, this free parameter is a remainder of O2- 
symmetry in the pure higgsino limit, and it does not enter into expressions for the masses 
(the last assertion is checked by direct calculation in the second order approximation). 

The structure of the neutralino chiral fields follows from the transformations 

0L = UXL, (pR = {(pLf = V*XR, 

XL = U-Vl = U+0i, Xr = iXif = U^0fl, (4.15) 

where U = U2 ■ U/ and (0)^ = {B,W^, H'^, H^). With the help of Eq. dUSj) for the 
non-chiral neutralino field X = Xl + Xi?; we get 

(p = {ReV - i75/mU)x, X = (ReV^ + 775 ImV'^)(j). (4.16) 

In the first order of mixing (see Appendix B) the structure of the neutralino fields is defined 
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by the following expressions 



Xi^ B + ^ + TT^hs^^f + 

V2 Ail + Ml - fi 



X2 



V^'M2+/i M2-/i 



)75i^2", (4.17) 



^3 - tJ^^^^ - I^^^ + + ^^^^ 

Mi+ fi M2 + fi V2 V 2 



X4 



Mi-/i M2-fi v/2 ' v/2 



Thus, the imaginary part of the transformations contains the factor corresponding to the 
redefinition x = ^TsX the minimal diagonalization procedure. It was checked by direct 
calculation up to the second order that the diagonalization of the real matrix with 
redefinition of the field with negative mass gives the same results when the free parameter 
is equal to zero (see Appendix B). 

It is known that in a wide class of SUSY scenarios the values of Mi 2 and/or fi are of the 
order of TeV and higher, so the coefficients in fl4.17p are of the order of 10^^ or less. Thus, 
the mixing terms give a contribution to the physical values ~ 1%, so the expressions (I4.17p 
can be used in practical calculations with a good accuracy. 



V. CONCLUSION 

It is known that the diagonalization of the neutralino mass form by the orthogonal real 
matrix leads to the neutralino mass spectrum with one negative mass. This has to be 
taken into account in calculation rules or by a redefinition of the field with negative mass. 
An alternative way is the diagonalization by a unitary complex matrix which leads to the 
mass spectrum with all positive masses. Formally, both the ways are equivalent, but the 
second one is more convenient, because it does not demand any modification of the standard 
calculation rules. 

In this work, we have considered the connection between the mass sign, the relative parity 
of the neutralino states and the structure of the boson-neutralino interaction. These features 
should be considered in the evaluation of the SI and SD contribution to neutralino-nucleon 
scattering. The suggested approach directly illustrates the existence of one free parameter, 
generated by the specific symmetry of the /i-term. When this parameter is equal to zero, 
both approaches give the same results. This was strictly shown in our work for the mass 
spectrum and states up to the second approximation. 

We suggest a simple and convenient way of diagonalization by a unitary matrix with 
the exponential parametrization. Having used this matrix, we get transparent perturba- 
tive formalization of the diagonalization procedure. This method gives simple expressions, 
illustrating the neutralino states structure and the form of the gauge boson-neutralino in- 
teraction. These expressions can be used in most of the SUSY scenarios with the accuracy 
~ 1% or higher. 
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Appendix A 

To explicitly show the appearance of the imaginary unit in the neutralino mass matrix and 
for consistency we give here the minimal part of the SUSY Lagrangian and briefly describe 
the transformation of the initial SUSY expressions to the ones in terms of four- dimensional 
fields. All definitions and calculations are in the notation of Refs. and j3l|. We consider 
the electro-weak part of the MSSM Lagrangian 

L = LG + LH + Lph. (5.1) 

In Eq.f l5.ip . the gauge term Lq has the standard form 

La = \{iW''W.)ee + mw'y.e + (W.^W'jee + {W^Wneg}, (5.2) 

where Wa and Wa are U{1) gauge superfields and Wj^ and are SU{2) gauge superfields. 
The Higgs term contains two chiral superfields with hypercharges Yi,2 = ±1, 

Lh = {H+ exp (^iGi - g2G2)Hi + 

+ H+exp {-giGi - g2G2)H2]eem- (5-3) 

The phenomenological part contains the so-called /i-term and gauge soft mass terms: 
Lph = ^i[{H^eH2)ee + {HteHt)ee\ - ^M^ibb + bb) - 

- ^M2{uJaUJa+U}aUJa), (5.4) 

where e = ir2. To define the notation of the components we also present the expressions for 
the gauge superfields Gi and G2 in Wess-Zumino gauge: 



Gi = Oa^e ■ + iOe ■ 9b - i99 ■ 9b + -99 ■ 99 ■ D^, 



G^ = 9(r''9 ■ Vr; + i99 ■ 9uj^ - i99 ■ 900" + ]^99 ■ 99 ■ D^, (5.5) 

W^ = -^DDD^G, Wa = -\DDD^G. 

The Higgs chiral superfields are 

Hi = hu + V29hi + 99 ■ Fi, H2 = hd + V29h2 + 99 ■ F2. (5.6) 
Thus, the particle content is 

Gi = (i?^,6),G^=(iy;,^"), 

H, = {hu,h),H2 = {hd,h2), (5.7) 
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where 

''« = (t°)-'" = (A'r)'"'=OI)'"^ = ('l)- 

In Eqs. fl5.7p and (15. 8p Bf^,WI^, hu, hd are boson fields and b,uj"-,hi,h2 are two-component 
fermions. 

Using the standard method from Eqs. fl5.1l) - fl5.4p we get Lagrangians in terms of two- 
component fermions. The gauge field Lagrangian is 

Lg = - \b^^'B,, + zba^^d.b + ^-Dl - \wrw;, 

+ lu'^a^id^UJa + g2eabc(^cWl) + ^/^s^^Sa, (5.9) 

where B^, = d^B, - d,B^ and W^, = d^W^ - d,W^ + g2e'''"'W^bW,c. The Higgs field 
Lagrangian is 

- ]^giB^hia^'hi + ]^giB^h2a^h2+ 

+ \g2W;h^a^T,h^ + \g2W;h2a^TM+ 

+ '-^{Kbhi - h,hK) - '-^iKbh2 - hbhd)- (5.10) 

- ^(/i+tu/ii - hiLohu) - ^{hjujh2 - h2tohd). 
The phenomenological Lagrangian is 

Lph = -^i{hieh2 + hieh2) -]^Mi{bb ^bb) - 

- ]^M2{uJai^2+i^^)- (5.11) 

In Eqs. fl5.9l) - fl5.11l) all fermion fields are two-component spinors. The transition to 
four-component Mayorana spinors in a Weyl basis is defined by the following relations 

Xk = (^''^ , 4>k = {b, uj", hi, /la); 

XkL = (^^Q^ , XkR = ^) , x1 = XR = Rx; (5.12) 

-5a--)-(JS)^4'^^-)-(o; 
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With fl512D we get 



ibaf'df.b = ^-B-i^d^B, bb + bb = BB, B = ^ 
iuo'^a^'d^Qa = ^W^V^M^" + Div{W), 

ha''h = HiYH^L, h2a^h2 = H2i^H2L] (5.13) 



- hibhu = h^BHiL - HilBK, 
h+bh2 - h2bhd = h+BH2L - H2LBhd; 
htujhi - IiiujK = KWHiL - HilWK, 
hjujh2 - h^ha = hjWH2L - H2LWhd; 
hieh2 + hieh2 = HiReH2L + HilsH2r. 

From Eqs. fl5.9p - fl5.lip with the help of Eqs. fl5.12p . fl5.13p we obtain Lagrangians in terms 
of four-component spinors. The gauge field Lagrangian is 

Lg = - Ib^'b,, + '-B^d.B + ^-Dl - \wrw;^ 

+ '-W'^Yid^Wa + 92eabcW'^W'^) + ^D^D2a. (5.14) 
The Higgs fermion field Lagrangian is 

Lh = '-HiY{df.H, + tgiB^H.L - ig2W^TaHiL) 
+ '-H2Y{d^H2 - tgiB^H2L - tg2W;TaH2L) 

+ '-^(h+BH,L - S,lBK) - '^(hjBH2L - H2LBha) (5.15) 

- '-^{h+WH,L - HilWK) - '-^{h+WH2L - H2LWhd). 

The phenomenological Lagrangian is 

LpH = -^M^BB - ^M2WaW^ - fi{H,ReH2L + HiLeH2R). (5.16) 
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Appendix B 



Here we consider a simple and easily formalized method of the complex mass form diag- 
onalization: 

(0^)^M,0i + h.c. = ixnfVjMjViXL + h.c. = niiXiXi- (5-17) 

In Eq. (15.171) M/ is a symmetric complex matrix, (p^ju^ are the chiral components of the 
initial Majorana spinor fields arising after intermediate diagonalization fl4.3p . and x ^^re the 
final Majorana fields (neutralino): 



lI\T 



{B,W\<Pi<l)i), X^ = {Xi,X2,X3,X4). (5.18) 



The intermediate states 03 and (pl are defined by Eq. (3.2). We suggest straightforward 
diagonalization of the form (I5.17P by the unitary matrix in the exponential parametrization. 
In the general case, the unitary matrix \J{n x n) has 2n^ — = parameters, where 
unitary conditions are taken into account. For n = 4 we have 16 input parameters, six 



angles and ten phases [12|. However, in the case of a symmetric mass matrix with real Mi^2 
and /i, this number of parameters is excessive, so we suggest a unitary matrix with six angle 
and six phase input parameters. These 12 parameters can be defined from 12 independent 
conditions following from the symmetric matrix diagonalization (Uf M/U/)it = 0, i > k (or 
i < k). It is convenient for the analysis to use U/ in the exponential form [32|: 



(5.19) 



In Eq. fl5.19l) 5i — 5q and 0i — 06 are the angle and phase parameters, respectively. The 
quantities 5^ and 0^ are the input parameters, while ap^r^^ak are some functions of the 
input parameters which follow from the unitary condition U"'"U = 1: 

ar = {l-6l- 61 - ^3^)^^ = - 6^ - rlY'\ 

as = (1 - 61 - rl - rlfl\ a, = {1 - rl - rl - rlfl\ 

airie^"^ + aa^ie^"^^ + 826 ^e"'^'''^-'''^'^ + 5^5r,e'^'*'^-^'^'^ = 0, 

airae^"^ + 5ir4e^("^+'^i) + a362e''^^ + d^dQc'^^-'-^"'^ = 0, (5.20) 

air3e^°3 + ^irgc^^^^+'^i) + (52r6e*('^^+°«) + 04536^'^^ = 0, 

^^^3ei(«3-"i) + asrse*"^ + (54r6e*('^*+°«) + 04556^'^^ = 0, 

^^^2gi(a2-ai) ^ a2r4e^°^ + 03^46^'^* + 656^6'^^'-^'^ = 0, 

r2r3e^("^-"^) + r^r^e'^'^'-'^*^ + a^r^e'""- + aid^e'^' = 0. 

The Anzats fl5.19p is convenient for approximate calculations in the case of a quasidiagonal 
mass matrix, for instance, M/ defined by Eq. (14.40 . In the case considered, the absolute 
values of the diagonal elements and the differences are much greater than the off-diagonal 
ones (the equality of the third and fourth diagonal elements fi in M/ is compensated by off- 
diagonal zero). The diagonalyzing matrix has a similar structure, i.e. the input parameters 
6k in Eq. (I5.20p are small, 6k ^ I and ~ 1. So, due to the smallness of parameters 6k ~ 
Mz/Ma, where a = 1, 2, one can easily solve the system of equations (I5.20p approximately. 
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The functions aa,rk and ak are determined from Eqs. fl5.20p . the input parameters 5k and 
(pk are defined by the conditions (UJMt-U/)^^ = 0, i > k. Hence, the diagonal elements 
(Uf M/U/)fcfc = nik give the masses in terms of known quantities. 

Finally, we have done the calculations up to the second order ~ M|/M^ (or M|//i^) 
inclusively and get the expressions for the elements of the diagonalizing matrix (15.191) (the 
hierarchy of and /i is arbitrary). Such an approximation is reasonable for calculations 
within a wide class of Split SUSY models. 

The input parameters are 



6ie 



5^e 



Ml sin 2ew M2 + /i sin 2/3 
2(M| - /i2) M1-M2 ' 



M2 + fi 
i ( 17111712 



636 



5 



+ 



M2 -/i 



2/i \Mi - fi M2- fi 



(5.21) 



The diagonal elements ap are 



ai 



a2 



0-3 



04 



777T 



+ 



rrin 



(Mi+/i)2 (Mi-/i) 



+ 



■777^ 



L(M2 + /i)2 (M2-/i)2j ' 



777T 



777n 



(Mi+/i)2 (M2+/i)^ 



+ 



777^ 



L(Mi-/i)2 (M2-/i)2 



(5.22) 



The off-diagonal elements are 



Tie 



r2e'' 



r4e 



104 



i«6 



+ 



77727774 



(Ml + /X) (M2 + /i) (Ml - /i) (M2 - /i) 

M| sin 2^vf(M2 + /isin2/?) 
2(M| -At2)(Mi -M2) ' 

«777i 7772 



Ml +/i ' 
~M2 + /i 

77717772 
Mf - 



Ml 



77737774 



7774 



Mo 



/i 



(5.23) 



+ 



M| 



11111712 



+ 



77737774 



2/i V^i - M2 - /i 



In Eqs. (15.211 - 15.231) we give the zero value to the free parameter Sg arising in the first 
order. It was checked in the first approximation that the existence of the free parameter Sq 
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leads to the phase redefinition of the fields xs and X4- Having applied Eqs. (15.211) - (15.231) . 
we obtain expressions for the neutralino masses (I4.14p : 

iVjMjVj)kk = rrik (5.24) 

which are positive. We have checked also the unitary condition U/U/ = 1. 
The structure of the neutralino chiral fields results from the transformations 

XL = U-Vl = U+0i, xr = iXLf = U^'^pR, (5.25) 

where U = Us ■ U/, Ua is defined by (EID and {(j)f = {B, W^, H^, H^). 

To illustrate the relation between the initial and physical fields we give transformations 
in the first order of mixing: 

imi m2 
Ml + /i Ml — II 

o irrij, 1714 

^y2 Mi+ fl Ml - fi 

i m4 1 i 

v/2 M2 + yu M2- n V2 V2 

„-o _ ^ ^ "^1 ^2 ^ _ 

V2 Ml + /i Ml - /i 

i m3 m4 1 i 

y/2 M2 + fi M2 - /i V2 V2 

Transformation of the i?-component can easily be found from the relation xr = {Xl)'~^ ■ 
Inverse transformations illustrate the neutralino structure: 

V2 Ml + /i Ml - /i 
v/2Wi+/i Mi-/i^ 

V 2 M2 + /U M2 - /i 
a/2 M2 + /i M2 - /i 



^mi , ^"^3 .,,3 I 1 0-0 , 1 0-0 

X4L ~ -77 -Bl + — -Wl - -^^IL + 



Mi-fi M2-/i ^ V2 72 



The transformation of the non-chiral Majorana field X = Xl + XR 

= (i?eU - i75/mU)x, X = {ReV^ + ^5 /mU^)0. (5.28) 
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In the first order of mixing from fl5.28p it follows that 



irrii 1712 

^ ^Xl- -Tl—, 75X3 - JJ X4, 

Ml + /i Ml — II 

^ ~ ^2 + , . . 75X3 + -TJ X4, 

M2 + /i M2 - /i 

i m3 m4 1 i 

+ ^ + JJ 75X2 + ^X3 - ^75X4, 

V2 M2 + /i M2- fx V2 V 2 

uo i , mi m2 

^2 ~ — 7^ ^fZ^ 17 75X1+ 

V2 Ml + /i Ml - fi 

i 1714 1 i 

+ ^ilTi—, JJ 75X2 + ^X3 + ^75X4- 

y/2 M2 + jj, M2- IJ, V2 V2 

The structure of non-chiral neutralino states is 

^ i , mi m2 s 

V2 Ml + /i Ml - /i 

+ v^Wi+/i Mi-zi^^^''^' 

V2 M2 + M2 - n 
Mi + fi M2 + fi V2 V 2 

Ml- n M2 - /i V2 V 2 

Thus, the imaginary part of the transformations contains the factor corresponding to the 
redefinition x = ^75X the intermediate diagonalization procedure. By direct calculation 
(up to the second order) it was checked that the diagonalization of the real matrix Nlj with 
redefinition of the final field with negative mass gives the same results as in our case with 
Sq = 0. Note that our formulae ( 15.261 - I5.30p coincide with the corresponding ones from 



Ref. [12| if we redefine the neutralino state with negative mass as x ^ '^75X- 

Now we generalize the method of diagonalization for the case of a mass matrix with 
complex parameters Mie*^^, M2e*'^2 and fie^^i^. Then we have to extend (15.191) introducing 
additional phase parameters according to O/j — >• a^e^^''. The functions ap^r^ and au are 
determined in terms of the input parameters 5k-, 4>k and by the unitary condition U^U = 1. 
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We represent them in the form 

ai = (1 - Si - Si - SlY'\ a, = {l-Sl- Sj - rlf/\ 
as = (1 - Si - rl - rlf'\ a, = (1 - rl - rl - rlf'\ 
a^^^ei(«i+ei) + a^^ie^^^i-^^) + S2Sie'^'''^-^*'^ + 
+ S^S^e^^^-'^'^ = 0, 

a^r2e'^"''+^') + 5ir4e*("*+'^^^ + a3S2e'^^''-^^^ + 

+ S^Sqc'^^''-'^''^ = 0, (5.31) 
+ a453e^('^3-^*) = 0, 

+ 55^66^^'^^"'^'^) = 0, 

^^^3e^(°3-ai) + a2r5e^("5+^^) + 54r6e^^'^*+"«)+ 
+ 04556^^'^^-^^) = 0, 

r2r3e^("^-"2) + r4r5e^("^-"^) + a3r6e^("«+^^)+ 
+ aiSee'^''''-^'^ = 0. 

The set of input parameters Sk, 4>k and .^^3 is determined in terms of the mass matrix elements 
utihzing the diagonahzation conditions (U^MU)^^ = 0, i ^ k and /m(U^MU)ij = 0. 

In the first approximation from the second condition we get ^1 = ■0i/2,^2 = '02/2 and 
^3 = ^4 = Prom the first condition in the same approximation we get 

5i = or 01 = Vi/2, Ml = M2; 



—mi 
M2 - 



[M2 + 2/xMi cos(V^i + V^^)]^/2^ 



tan(</>2 - y) = cot(^— ); 



-m2 
Ml-n 



[Ml + + 2/xMi cos(0.i + 0.^)] 



tan(03 - y) = tan(^— ); (5.32) 

X "^3 r,^2 , .2 



[M2^ + _ 2^M2 COS(02 + 0M)]^^^ 



tan(04 - y) = cot(^^); 

, [M^ + ^ 2/XM2 cos(0;2 + 0^^)! 



tan(0, - y) = tan(^— ); 

^6 - or 06 = (02 + 0'm)/4- 
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All expressions for the input phases can be represented as explicit functions, for example 

02 = Y - arctan( ^^ cot( ^ )). (5.33) 

From Eqs. (15.321) one can see that taking account of the complex degrees of freedom in the 
mass matrix complicates the calculations. However, the perturbation scheme of the method 
is retained and can easily be formalized. It has to be noted also that we use the same number 
of phases as in the general method [12I], and the question of optimization of this number for 
a special kind of mass matrix is practically important. 
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